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I.  INTRODUCTION 


In  principle,  the  concept  for  determining  the  optimum  location  of  an 
electromagnetic  sensor  on  or  near  an  airborne  platform  such  that  the  sensor 
responds  to  the  incident  field  and  not  to  the  total  field  is  simple.  All  one 
needs  to  do  is  to  determine  a  position  for  the  sensor  where  the  scattered 
component  is  minimum  for  all  possible  excitations  and  frequencies  and  locate 
the  sensor  at  that  position.  If  the  sensor  is  polarization  sensitive,  only 
the  responding  components  of  the  scattered  field  need  to  be  determined.  The 
problem  could  be  approached  analytically,  numerically,  experimental ly,  or  any 
combination  thereof.  Here  we  present  results  of  the  scale  model  experimental 
study  carried  out  in  the  University  of  Michigan  Radiation  Laboratory  near 
field  measurement  facility  using  1:72  scale  models  of  F-106  aircraft. 

In  the  anechoic  chamber  measurements  of  signal  picked  up  by  a  miniature 
sensor  that  represents  the  boom  B-dot  sensor  were  made  with  and  without  the 
F-106  model  present  in  both  the  ground  plane  (Sec.  3)  and  the  free  space 
configurations  (Sec.  4).  In  total,  17  different  excitation  situations  were 
simulated.  For  each  excitation  nine  points  in  the  frontal  region  of  the  model 
were  probed  over  118-4400  MHz  corresponding  to  1.64-61.3  MHz  at  full  scale. 
Since  the  main  aircraft  resonances  occur  at  7.4  MHz  and  13.7  MHz  for  the 
fuselage  and  the  wings,  respectively,  the  resonances  fall  well  within  the 
frequency  range  covered  in  the  measurements.  For  each  excitation  and  measure¬ 
ment  location  (153  in  total)  an  error  due  to  scattering  from  the  model  is 
determined.  In  general,  the  results  agree  with  expectations:  the  further  the 
sensor  is  away  from  the  aircraft  the  smaller  is  the  error  and  the  error  is 
least  when  the  sensor  is  located  in  the  plane  of  the  wings  of  the  aircraft. 
The  important  contribution  of  this  study  is  the  determination  of  actual  errors 
that  can  be  anticipated  for  the  given  sensor  location  and  excitation.  In 
general,  these  errors  are  small,  mainly  due  to  the  almost  ideal  geometry  of 
the  F-106  aircraft  that  has  minimal  scattering  in  the  forward  region.  The 
assessment  and  discussion  of  errors  is  presented  in  Section  5. 

Other  studies  performed  under  this  study  include  (a)  SEM  studies  of  F-106 
aircraft,  and  (b)  lightning  simulation  studies  in  the  frequency  domain.  These 
are  presented  in  the  appendices. 
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II.  FACILITIES,  INITIAL  STUDIES  AND  OPTIMIZATION  CRITERIA 


1.  MEASUREMENT  FACILITY  AND  INSTRUMENTATION 

The  University  of  Michigan  field  measurement  facility  (Ref.  1),  that  has 
been  used  almost  exclusively  for  measurements  of  aircraft  exterior  field 
coupling  to  the  incident  electromagnetic  field  (Refs.  2,  3),  is  "tailor  made" 
for  the  present  study.  A  block  diagram  of  this  facility  is  shown  in  Figure 
1.  The  system  is  a  CW  one  in  which  the  frequency  is  swept  (stepped)  over  a 
wide  range.  A  key  part  of  the  facility  is  a  tapered  anechoic  chamber  approxi¬ 
mately  50  feet  in  length  with  the  rectangular  test  region  18  feet  wide  and  12 
feet  high.  Its  rear  wall  is  covered  with  72-inch  high-performance  pyramidal 
absorber,  with  18-inch  material  used  on  the  side  walls,  ceiling  and  most  of 
the  floor.  The  material  in  the  tapered  section  (or  throat)  is  two-inch  hair- 
flex  absorber.  The  chamber  can  be  thought  of  as  a  lossy-wall  horn  antenna 
terminated  by  the  rear  wall.  The  signal  is  launched  from  a  single  exponen¬ 
tially  tapered  broadband  antenna  located  at  the  apex  of  the  chamber.  The 
antenna  is  fixed,  and  since  the  radiated  signal  is  horizontally  polarized,  the 
pseudoplane  wave  in  the  center  ('quiet  zone')  portion  of  the  test  region  is 
also  horizontally  polarized. 

The  instrumentation  is  centered  around  a  Hewlett-Packard  84108  network 
analyzer,  and  is  computer  controlled.  An  HP9830A  computer  controls  the  fre¬ 
quency  to  be  generated,  switches  in  the  appropriate  power  amplifiers  and  low- 
pass  filters,  and  reads  and  stores  the  amplitude  and  phase  of  the  signal 
picked  up  by  the  sensor.  During  a  run,  the  frequencies  are  typically  stepped 
from  118  to  4400  MHz.  Because  of  the  limited  memory  size  of  the  computer, 
this  frequency  range  is  recorded  in  four  bands:  118  to  550  MHz  (in  4.8  MHz 
steps),  550  to  1100  MHz  (in  4.8  MHz  steps),  1100  to  2000  MHz  (in  9.6  MHz 


1.  EMP  Interaction:  Principles,  Techniques  and  Reference  Data,  AFWL-TR-80- 
402,  pp.  267-2^5,  December  1^80,  EMP  Interaction  Note  2-1. 

2.  Liepa,  V.V.,  Y.  Hyun,  "Free  Space  Scale  Model  Measurements  of  the  F-14A," 
University  of  Michigan  Radiation  Laboratory  Report  No.  017463-4-T;  Inter¬ 
action  Application  Memo  38,  March  1982. 

3.  Liepa,  V.V.,  "Free  Space  Mode  FB-111A  Scale  Model  Measurements,"  Univer¬ 
sity  of  Michigan  Radiation  Laboratory  Report  No.  017463-5-T;  Interaction 
Application  Memo  35,  July  1982. 
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steps),  and  2000  to  4400  MHz  (in  16.8  MHz  steps).  The  data  from  each  band  are 
stored  by  the  HP9830A  computer  on  a  cassette  for  later  transfer  to  an  HP9845B 
computer  which  processes  and  plots  the  data.  If  substantial  processing  or 
computation  is  involved,  or  if  a  need  exists  to  write  the  data  on  standard 
computer  tape,  the  data  are  transmitted  to  the  central  University  of  Michigan 
AMDAHL/5860  (IBM  compatible)  computer. 

The  signal  measured  by  the  network  analyzer  as  picked  up  by  a  miniature 
probe  or  a  sensor  is  a  function  not  only  of  the  measured  field  but  of  the 
entire  facility,  including  the  probe,  chamber,  antenna,  amplifiers  and 
cabling,  and  it  would  be  virtually  impossible  to  separate  and  to  correct  for 
the  contributions  of  each.  The  approach  taken  in  this  study  is  to  apply  an 
appropriate  calibration  or  normalization  whereby  the  response  of  the  facility 
and  the  instrumentation  are  eliminated.  In  short,  a  measurement  is  made  with 
an  airplane  model  present  (test  measurement)  and  then  repeated  with  the  air¬ 
plane  model  removed  (calibration).  Taking  the  ratio  of  the  measurements 
produces  the  required  data  and  cancels  out  the  facility  responses  and  measure- 

Jr 

ment  instrumentation  responses. 

Using  the  fundamental  properties  of  electromagnetic  scattering  by  simple 
symmetric  bodies,  Gi  ri  and  Baum  (Ref.  4)  have  provided  guidelines  for  place¬ 
ment  of  the  sensors  on  or  near  the  aircraft,  but  due  to  their  nonexact  analy¬ 
sis,  their  error  estimates  are  qualitative  for  most  of  the  frequency  ranges  of 
interest.  Using  their  guidelines  and  the  physical  constraints  imposed  by  the 
safety  and  performance  consideration  of  the  F-106  aircraft,  the  B-dot  sensor 
location  was  limited  to  the  frontal  regions  of  the  aircrft  with  the  sensor 
axis  parallel  to  the  fuselage  and  in  the  plane  of  symmetry  of  the  aircraft. 

Our  task  was  to  find  the  optimum  position  for  the  sensor,  limited  to  the 
frontal  region  of  the  aircraft.  Figure  2  shows  a  side  view  of  the  aircraft 
with  the  nine  measurement  locations  indicated  by  x,  where  the  measurements 
were  made  for  specified  illuminations  with  and  without  the  model  present. 
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Gi ri ,  O.V.  and  C.E.  Baum,  "Airborne  Platform  for  Measurement  of  Transient 
of  Broadband  CW  Electromagnetic  Fields,"  Sensor  and  Simulation  Note  284, 
May  1984. 


These  locations  were  chosen  to  cover  the  area  within  which  the  sensor  place¬ 
ment  is  feasible  and  have  been  (theoretically)  investigated  by  Giri  and  Sands 
(Ref.  5). 


2.  INITIAL  STUDIES 

In  order  to  better  understand  the  results  of  nose  area  measurements  and 
interpret  the  various  oscillations  in  the  frequency  response,  surface  current 
measurements  were  first  made  to  determine  the  fuselage  and  the  wing  resonances 
of  the  F-106A  using  small  current  loops  (Refs.  2,3).  Figure  3  shows  the 
normalized  axial  current  density  (J_  =  n  x  ]1,  A/m)  on  top  of  the  fuselage  for 
the  case  of  top  incidence.  E-parallel  to  the  fuselage.  The  location  is  at 
midpoint,  measured  horizontally  halfway  between  the  bulkhead  and  the  top  of 
the  vertical  fin.  The  (top)  curve  shows  the  fuselage  resonates  at  7.2  MHz  and 
reaches  a  value  of  about  10.5  at  the  peak  (relative  to  the  incident  field). 
Note  the  model  used  was  without  the  Pitot  boom  and  the  radome.  The  lower 
curve  shows  the  corresponding  phase,  with  reference  plane  at  the  point  of  the 
measurement.  The  physical  optics  (P.0.)  estimates  are  indicated  by  the  dashed 
lines,  which  are  twice  the  incident  (magnetic)  field  for  the  current  and  zero 
degrees  for  the  phase. 

Figure  4  shows  the  wing  current  measured  on  the  bottom  of  the  fuselage 
along  the  line  of  the  flap  joins,  for  bottom  incidence.  E-perpendicular  to  the 
fuselage.  In  the  model  the  flaps  were  rigid,  without  a  gap  between  the  wing 
and  the  flap.  Note  the  wing  resonance  occurs  at  13.7  MHz,  with  the  maximum 
current  reaching  about  3.5. 

The  basic  dimensions  of  the  F-106A  are  shown  in  Figure  5.  Using  the 
fuselage  length  of  17.67  m  (bulkhead  to  tail)  and  the  resonant  frequency  of 

7.4  MHz,  the  fuselage  resonance  length  is  0.436  wavelength,  which  is  as 

expected  for  a  relatively  fat  body.  The  wing  resonance,  however,  occurs  when 
the  wing  span  is  0.533  wavelength  (wing  span  11.67  m,  resonant  frequency 
13.7  MHz);  this  contradicts  the  basic  past  observations  that  most  of  the  main 
resonances  occur  when  the  body  dimensions  are  0.4  to  1.47  wavelengths. 

Perhaps  in  this  case  the  resonance  observed  is  not  solely  due  to  the  mode 
supported  by  the  wings  but  characterized  by  the  wing  and  the  vertical  fin 

characteristic  dimensions. 
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Giri,  D.V.,  and  S.H.  Sands,  Design  of  Incident  Field  3-0ot  Sensors  for 
the  Nose  8oom  of  the  NASA  F- 1068  Aircraft,  AFWL-TR-^3- 14 1,  April  1984. 


FREQUENCY  ( MHz ) 


Surface  current  for  F-106A  on  top  of  the  fuselage,  top 
incidence  E-parallel  to  fuselage. 
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3.  DATA  CONDITIONING 


The  data  plotted  in  Figures  3  and  4  are  as  measured  without  any  smoothing 
or  filtering.  As  seen,  it  is  somewhat  noisy  at  the  lower  end,  and  shows  some 
drop  outs  and  small  oscillations  that  definitely  are  not  a  part  of  aircraft 
response.  Hence,  some  smoothing  is  justifiable  and  when  done  properly  will 

improve  the  accuracy  as  well  as  the  appearance  of  the  data.  All  the  data  used 

for  optimization  of  B-dot  sensor  was  filtered  using  program  FILTR5*  that  is 
based  on  the  convolution  principle  and  convolves  a  sine  (=(sinx)/x)  function 
with  the  frequency  data.  The  convolution  of  the  sine  function  with  frequency 
is  equivalent  to  gating  or  multiplying  the  equivalent  time  domain  signal  by  a 
window  function.  For  filtering  the  data  (Sec.  3  and  4),  30  cycles  were 

selected  as  the  cut-off  frequency  for  the  oscillations  on  the  frequency 
data.  This  is  equivalent  to  a  range  gate  of  about  130  m  or  six  aircraft 
lengths.  In  relation  to  model  measurements  (1/71.8  scale  modeling)  this 

corresponds  to  a  distance  of  1.8  m  which  is  also  the  distance  from  the  model 
to  the  ceiling  or  to  the  floor  in  the  chamber.  The  model  to  side  wall  dis¬ 
tance  is  about  2.2  m  and  hence  filtering  also  removes  these  interactions  if 
such  exists. 

On  the  other  hand,  if  a  model  shows  a  high  Q  response,  filtering  can  also 
decrease  the  resonant  peak  and  create  erroneous  data.  Note,  a  resonance  can 
be  treated  as  a  signal  bouncing  back  and  forth  on  the  body,  and  if  this  bounc¬ 
ing  continues  longer  than  the  time  window,  a  part  of  the  response  would  be  cut 
off  resulting  in  a  lower  peak.  It  is  therefore  essential  that  filtering  be 
tested  on  the  worst  case  data  -  high  noise  and  high  resonance  peaks  (or  high  Q 
response)  to  determine  the  acceptable  filtering  parameters. 

Figure  6  shows  the  worst  case  data  that  was  obtained  in  this  study.  The 
measurement  has  rather  high  (chamber)  interference,  dropout  noise,  and  the 
highest  model  interaction  (i.e.,  deviation  from  unity )  observed  of  all  the 
situations  measured.  The  corresponding  filtered  plot  is  shown  in  Figure  7. 
Note,  the  noise,  the  spikes,  and  the  fast  wiggles  have  disappeared,  but  the 
major  features  that  arise  from  the  model  interaction  have  remained  intact. 


In-house  program  developed  for  HP9845T  to  process  the  data. 


4.  ERROR  EVALUATION  CRITERION 

Depending  on  direction  of  Illumination,  polarization,  and  the  location  of 
the  sensor,  the  data  recorded  varies  from  the  worst  case  shown  in  Figure  7  to 
almost  a  straight  line  (no  interaction).  In  total  there  were  156  curves 
generated  and  in  order  to  categorize  or  evaluate  these  curves,  error  values  Ej 
and  E2  are  defined  giving  a  pair  of  numbers  for  each  data  curve.  These  error 
definitions  follow  from  average  and  r ms  function  definitions  and  are  defined 
(in  percent)  as 


where 

N  *  the  total  number  of  data  points  (118-4400  MHz,  N  -  450  ) 
fp  *  the  measurement  frequency 
aF  »  the  total  frequency  interval 
Af n  »  fn  -  fn_1.  sampling  interval 

v(fp)  *  the  signal  measured  with  the  model  present 
v^fp)  s  the  signal  measured  without  the  model 

vo^n^  3  si9nal  measured  without  the  model,  but  the  sensor  rotated 

for  maximum  signal  and 

K  *  phase  correction  constant 

Note  that  in  a  case  such  as  the  top  incidence,  H-parallel  to  the  fuselage 
v0(fn)  3  vl^n^»  whereas  f°r  t0P  incidence.  E-parallel  to  the  fuselage, 
v0(fp)  f  Vi(fp).  In  fact,  vj(fp)  =  0,  and  is  only  non-zero  due  to  the 

reflection  of  the  signal  from  the  walls  of  the  anechoic  chamber,  the  support 
pedestal,  and  other  imperfections.  The  phase  term  K  was  added  to  allow 

correction  for  the  linear  change  in  the  phase  between  the  model  and  without 
the  model  data  that  results  from  dimensional  changes  in  the  chamber,  cabling 
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and  even  equipment  due  to  changes  in  temperature,  humidity,  etc.  To  minimize 
such  drift,  the  measurement  sequence  consisting  of  one  illumination  and  nine 
measurement  points  is  completed  in  the  same  day,  and  even  then  electrical 
changes  in  the  system  occurred  which  resulted  in  phase  drifts  on  the  order  of 
10  to  20  degrees.  These  can  be  related  to  a  few  millimeter  changes  in  the 
chamber  or  cable  dimensions.  In  the  time  domain,  a  linear  phase  shift  corre¬ 
sponds  to  the  time  delay  (or  advance)  with  no  change  in  pulse  shape,  and  hence 
one  is  justified  in  adjusting  the  linear  phase  term  K  to  minimize  the  error 
quantities  Ej  and  ^  Examination  of  the  phase  plot  in  Figure  7  shows  a 
slight  rise  in  phase  with  the  frequency  indicating  drift  in  the  equipment.  If 
not  corrected  for  this  slope,  E^  *  18  percent  and  when  corrected,  E1  =  11 
percent. 

For  processing  the  data  to  evaluate  E^  and  E2  and  optimize  the  phase  K,  a 
program  SUBTR5*,  based  on  Equations  1  and  2,  was  written  and  used. 


*  In-house  program  developed  for  HP9845T  to  process  the  data. 


III.  GROUND  PLANE  MEASUREMENTS 


When  a  model  to  be  measured  is  symmetric  about  a  plane  and  the  excitation 
is  such  that  incident  electric  field  vector  is  perpendicular  to  this  plane, 
measurements  can  be  made  using  one-half  the  model  mounted  on  a  large  metal 
plate.  The  plate  should  be  large  enough  to  avoid  errors  caused  by  edge 
diffraction  and  the  model -edge  interactions.  The  main  advantage  of  using  an 
image  or  ground  plane  technique  is  to  avoid  exposing  probe  leads  to  the  inci¬ 
dent  fields  to  interfere  with  the  measurements.  In  addition,  the  ground  plane 
has  the  tendency  to  improve  the  field  distribution  by  enforcing  that  the 
electric  field  be  perpendicular  to  the  metal  plate.  On  the  other  hand,  the 
disadvantage  of  using  a  ground  plane  is  that  it  limits  the  measurements  to 
symmetric  bodies  and  limits  the  illumination  to  the  case  of  incident  electric 
vector  perpendicular  to  the  ground  plane  if  the  free  space  plane  wave  inci¬ 
dence  is  to  be  simulated. 

Because  image  plane  measurements  do  provide,  in  general,  more  accurate 
results,  measurements  that  were  appropriate  to  such  configuration  were  per¬ 
formed  using  image  plane,  and  the  remaining  ones  were  performed  using  the  free 
space  technique  (see  Sec.  IV). 

1.  GROUND  PLANE  DESIGN 

As  mentioned  in  Section  II  the  anechoic  chamber  simulates  a  free  space 
environment  and,  since  it  would  require  considerable  modification  to  create  a 
ground  plane  chamber,  it  was  necessary  to  construct  a  ground  plane  that  could 
be  inserted  into  the  test  region  of  the  present  chamber  as  indicated  in  Figure 
1.  In  our  chamber  the  incident  E-field  propagation  direction  is  fixed  and  the 
E-field  is  horizontally  polarized.  It  is  therefore  necessary  to  have  the 
ground  plane  vertical,  with  the  capability  of  rotating  the  model  relative  to 
the  incident  field. 

For  the  reasons  given  below,  we  chose  a  circular  aluminum  ground  plane 
with  a  resistive  sheet  extension  as  shown  in  Figure  8.  The  aluminum  disk  is 
1.53  m  (5  ft)  in  diameter  and  is  supported  vertically  by  a  wooden  stand  at  a 
height  midway  between  the  floor  and  ceiling  of  the  chamber,  aligned  such  that 
the  incident  field  is  almost  edge-on  to  the  disk.  In  practice,  because  of  the 
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small  stray  (scattered)  signals  that  are  always  present  in  a  measurement 
situation,  the  ground  plane  was  positioned  5  deg  (counterclockwise  about  a 
vertical  axis)  to  insure  a  uniform  excitation  over  the  working  area  of  the 
disk. 

It  was  the  original  intent  to  permanently  attach  the  model  to  the  ground 
plane  and  to  produce  the  different  excitations  by  rotating  the  aluminum  disk 
about  a  horizontal  axis.  This  proved  to  be  impractical,  partly  because  of 
slight  warping  and  irregularities  in  the  ground  plane,  but  more  because  of  the 
sheer  size  and  weight  of  the  structure.  It  was  therefore  decided  to  keep  the 
disk  stationary  and  permit  the  rotation  of  the  model  in  discrete  steps. 

Although  the  disk  was  not  rotated  as  originally  planned,  it  is  believed 
that  the  circular  shape  is  still  desirable.  With  any  finite  ground  plane 
there  is  the  possibility  of  edge  waves  which  can  perturb  the  uniform  field 
characteristics  of  an  infinite  structure  and,  when  the  model  is  introduced,  an 
edge-model  interaction  can  occur  which  would  not  be  eliminated  by  the  field 
normalization  employed.  By  virtue  of  the  circular  geometry,  whatever  interac¬ 
tion  does  occur  should,  to  a  first  order,  be  independent  of  the  model  rota- 
ti  on. 

However,  for  a  given  illumination  a  situation  can  exist  where  a  traveling 
wave  can  be  excited  on  the  ground  plane.  This  wave  travels  in  the  direction 
of  the  incident  wave  (from  right  to  left  in  Fig.  8),  and  can  reflect  from  the 
left  edge  and  cause  disturbance  in  the  field.  To  minimize  such  reflections  a 
resistive  sheet  structure  was  built  and  attached  to  the  rear  edge  as  shown  in 
Figure  8.  The  sheet  was  curved  towards  the  back  to  hide  the  rear  edge  and  was 
constructed  by  applying  layers  of  resistive  paint  to  art  paper  supported  by  a 
wooden  frame.  By  varying  the  paint  thickness  the  resistivity  could  be  con¬ 
trolled,  and  the  final  structure  had  a  resistivity  which  started  at  about 

2  2 
10  2  at  the  metal  edge  and  increased  uniformly  to  about  1000  2  at  the  outer 

(left-hand)  edge. 

2.  GROUND  PLANE  MODEL 

The  ground  plane  model  was  constructed  from  a  1/72  scale  plastic  model 
(Hasegawa,  No.  1054)  which,  after  assembly,  was  cut  in  half  and  attached  to  a 
thin  aluminum  sheet  approximately  0.07  cm  thick.  The  model  was  then  sprayed 


with  silver  paint  to  simulate  the  metallic  characteristics  of  the  actual 
aircraft.  The  use  of  metal  plate  insured  ruggedness  which  was  necessary  for 
frequent  repetitive  use  of  the  model  on  the  ground  plane. 

To  mount  the  model  on  the  ground  plane,  aluminum  tape  was  used  to  tape 
the  plate  on  which  the  model  was  mounted.  This  tape  also  assured  a  good 
electrical  continuity  between  the  plate  and  the  ground  plane. 

Figure  9  illustrates  the  model  and  gives  its  dimensions.  The  actual 
scale  factor  of  the  model  is  1:71.8  and  is  based  on  the  ratio  of  the  wingspan 
dimensions  which  are  the  critical  ones  for  E-perpendi cular  to  the  fuselage 
polarization.  Figure  10  shows  photographs  of  the  model. 

3.  MEASUREMENT  PROCEDURE 

The  central  portion  of  Figure  8  depicts  the  model  and  the  probe  mounting 
on  the  ground  plane.  The  probe,  an  MGL-8,  is  mounted  at  the  center  of  the 
circular  ground  plane,  and  the  model  with  its  mounting  plate  is  repositioned 
on  the  ground  plane  to  attain  the  nine  probe  positions  relative  to  the  air¬ 
craft  as  given  in  Figure  2. 

The  measurements  were  made  for  incidence  of  (a)  Nose-on  (e  *  180  deg), 
(b)  Top  (9  =  -10  deg),  (c)  Tail-on  (e  =  0  deg),  (d)  Below  Tail-on  (9  =  -10 
deg)  and  (e)  Above  Tail-on  (9  =  10  deg).  For  cases  (a),  (b)  and  (c)  as  the 
model  was  rotated  so  was  the  probe  to  measure  the  magnetic  field  component 
parallel  to  the  fuselage.  The  cases  (d)  and  (e)  were  included  in  the  study  to 
verify  if  a  slight  tilt  of  the  sensor  axis  can  reduce  further  the  scattering 
errors  from  the  wings.  For  convenience  in  performing  the  measurements  the 
model  was  tilted  rather  than  the  probe  and  is  justifiable  due  to  the  broad 
(dipole)  pattern  of  the  sensor. 

In  addition  to  the  movement  of  the  model  relative  to  the  probe  (pos  1 
through  9)  and  the  rotation  of  the  model  with  the  sensor  to  obtain  the  inci¬ 
dences  (a)  through  (e)  given  above,  the  ground  plane  was  also  rotated  about 
its  veritical  axis  through  <p  =  5,  20,  45  and  60  deg.  The  <j>  =  5  degrees 
measurements  correspond  to  the  free  space  plane  wave  incidence  case,  and  the 
other  cases  could  imply  a  dual  plane  wave  incidence  exciting  only  the  anti¬ 
symmetric  current  components  on  the  model.  Figure  11  defines  the  9  and 
1)  variables. 


Photographs  of  the  ground  plane  model.  Top:  model  mounted 
on  a  thin  aluminum  plate.  Bottom:  the  same  mounted  on  the 
ground  plane  for  nose-on  incidence  H-vertical.  Note  the 
MGL-8  sensor  near  the  bulkhead. 
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As  a  test  of  the  facility  and  the  MGL-8  sensor,  the  sensor  was  first 
mounted  at  the  center  of  the  ground  plane  and  the  signal  recorded  as  a  func¬ 
tion  of  frequency  for  angles  of  0  from  0  to  360  deg  in  increments  of  15  deg. 
The  amplitude  curves  versus  0  exhibited  sinusoidal  variation,  but  there  was  a 
slight  shift  in  the  curve  implying  that  the  incident  signal  is  not  exactly 
propagating  horizontally,  but  on  the  average  arrived  along  the  line  0  =  -2 
deg.  This  could  be  a  result  of  antenna  placement  below  the  centerline  or  a 
difference  in  the  absorber  materials  between  the  top  and  the  bottom  in  the 
tapered  feed  region  of  the  chamber.  As  the  result  of  this  analysis,  the  e  =  0 
deg  specification  of  the  model  rotation  is  not  exactly  horizontal  but  tilted  2 
deg  clockwise. 

Consider  next  the  procedure  for  measuring  the  case  of  nose-on 
incidence  (0  =  180  deg  and  <j>  =  5  deg),  the  situation  for  which  the  layout  of 
Figure  8  is  applicable.  The  loop  axis  for  the  sensor  is  horizontal,  but  the 
incident  magnetic  vector  is  vertical,  a  case  where  there  should  be  no  signal 
(ideally)  received  by  the  sensor.  Next  refer  to  Equations  1  and  2  and  recall 
the  measurements  needed  to  compute  the  errors  Ej  and  E£.  These  are: 

v(fn)  =  sensor  output  voltage  with  the  model  present 

vi(f„)  =  sensor  output  voltage  without  the  model  (invisible  aircraft) 

vQ(fn)  3  sensor  output  voltage  with  the  sensor  rotated  to  measure  the 
maximum  magnetic  field  (without  aircraft) 

Also  refer  to  Figure  2  showing  the  placement  of  the  model  relative  to  the 
sensor.  Recall  that  the  sensor  is  fixed  and  the  model  is  moved  to  obtain  the 
needed  positioning.  For  the  ground  plane  measurements  the  numbers  2.5,  3.7 
and  5  cm  refer  to  the  distance  from  the  center  of  the  sensor  to  the  bulkhead 
of  the  model,  and  the  numbers  +0.5  and  -0.5  cm  refer  to  the  position  of  the 
model  relative  to  the  sensor.  Thus,  the  set  (5.0,  +0.5)  corresponds  to  the 
sensor  position  9  (Fig.  2).  (For  the  free  space  measurements  presented  in 
Section  4  the  convention  was  changed  whereby  the  numbers  +0.5  and  -0.5  cm 
refer  to  the  position  of  the  sensor  relative  to  the  model  and  hence  the  set 
[5.0,  +0.5]  would  correspond  to  the  sensor  pos  3.) 

In  carrying  out  the  measurement  sequence  the  voltages  v(fn)  are  first 
recorded  for  four  or  five  sensor  positions,  then  Vj(fn)  for  the  case  without 
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the  model,  and  then  followed  by  v(fn)  measurements  for  the  remaining  posi¬ 
tions.  If  the  sensor  needs  to  rotate  90  degrees  to  measure  the  incident  field 
v0(fn),  such  measurement  is  made  last.  To  minimize  the  measurements  errors 
due  to  equipment  drift  and  chamber  changes  due  to  the  temperature,  humidity, 
etc.,  the  measurements  sequence  is  completed  in  the  same  day. 

The  data  are  processed  using  program  SUBTR5  based  on  Equations  1  and  2. 
For  given  position,  say  position  1,  minimum  E^  is  computed  by  varying  the  K. 
The  search  was  done  manually,  and  typically  took  on  the  order  of  five  educated 
tries  to  arrive  at  min  with  K  increments  of  0.01  which  correspond  to  a  mere 
0.6  deg  ambiguity  at  the  maximum  measurement  frequency. 


4.  GROUND  PLANE  DATA 


The  results  are  presented  for  11  orientations  or  illuminations  of  the 
model.  For  each  situation  there  were  at  least  10  data  runs  generated,  9  for 
the  different  probe  positions  and  at  least  one  for  the  reference,  without  the 
model  present.  For  each  situation  studied,  a  summary  sheet  is  presented 
giving  Ej^  and  E2  error  values,  data  file  names,  reference  (without  airplane) 
file  names,  and  other  notes  that  may  be  essential  for  retrieval  of  the  data 
from  storage  tapes  for  further  analysis  analysis.  Each  summary  sheet  is 
followed  by  two  response  curves  for  sensor  positions  4  and  6,  respectively. 
These  give  a  descriptive  indication  of  the  behavior  of  the  received  signal  as 
a  function  of  frequency  for  the  particular  excitation.  Curves  for  the  other 
sensor  positions  were  not  included,  but  data  can  be  provided  to  subsequent 
users,  if  needed,  by  the  authors. 


Table  1  summarizes  the  data  presented  for  the  ground  plane  measurements. 


TABLE  1.  SUMMARY  OF  GROUND  PLANE  MEASUREMENTS 
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Set  3.3(b):  Normalized  field  at  position  4. 
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Set  3.5(b):  Normalized  field  at  position  4. 
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Set  3.6(c):  Normalized  field  at  position  6. 
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Set  3.7(c):  Normalized  field  at  position  6. 
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Set  3.10(a):  E^  and  E^  for  top  incidence,  H-parallel  to  fuselage,  i>  =  45  deg. 
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Set  3.10(b):  Normalized  field  at  position  4. 
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Set  3.11(b):  Normalized  field  at  position  4. 
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Set  3.11(c):  Normalized  field  at  position  6. 
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IV.  FREE  SPACE  MEASUREMENTS 


Not  all  of  the  incident  field  excitations  that  are  of  interest  can  be 
simulated  in  the  ground  plane  configuration,  which  basically  limits  the 
measurements  to  the  case  of  E-perpendicular  to  the  fuselage  or  antisymmetric 
excitations.  With  free  space  measurements  the  limitations  are  not  as  severe, 
but  there  is  a  possibility  that  an  external  sensor  and  its  metallic  leads  may 
cause  undue  interference  with  the  incident  field  and  the  model.  By  careful 
alignment  of  the  lead  normal  to  the  incident  electric  field  and  keeping  this 
lead  extending  away  from  the  model,  the  interaction  can  be  kept  down  to  a 
negligible  level. 

1.  MODEL,  SENSORS,  AND  MEASUREMENTS 

To  perform  the  free  space  measurements  the  ground  plane  was  removed  from 
the  chamber  and  the  Styrofoam  support  pedestal  which  is  used  to  support  the 
free  space  models  was  reinserted.  The  instrumentation  and  data  acquisition 
procedures  remained  essentially  the  same. 

The  same  size  (1:72  scale)  and  make  (Hasegawa,  No.  1054)  plastic  model 
was  built  and  prepared  as  discussed  in  Section  3.2,  except  it  was  not  cut  in 
half.  Intentionally,  the  fuselage  length  and  the  wingspan  were  adjusted  to 
dimensions  of  1/71.8  scale  the  same  as  for  the  ground  plane  case.  The  model 
is  that  of  the  F-106A  and  there  was  no  attempt  to  modify  it  into  an  F-106B 
version,  since  the  difference  in  the  two  versions  is  in  the  canopies  and  has 
little  effect  on  scattering  by  the  forward  region  where  the  fields  are  being 
studied. 

Two  free-space  probes  or  sensors  were  used  for  these  measurements:  (a) 
one  MGL-9(R)  and  (b)  our  own  make--a  3  mm  dia  shielded  loop  probe  made  of 
0.020-in-dia  coaxial.  For  measurements  either  probe  was  suspended  down  from 
the  center  into  the  chamber  such  that  the  coaxial  lead  was  perpendicular  to 
the  incident  electric  field.  The  MGF-9(R)  sensor  has  dual  coaxial  leads 
extending  radially,  and  to  adapt  it  to  a  situation  that  would  measure  the 
incident  H-field  and  at  the  same  time  have  the  lead  perpendicular  to  the  E- 
field  required  bending  the  leads  as  shown  in  Figure  12.  (An  axial  lead  ver¬ 
sion  of  the  sensor  MGL-9(A)  would  be  directly  applicable  but  such  was  not 
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Figure  12.  Modified  MGL-9(R).  The  coaxial  cables  were  bent  4.75  in 
from  the  loop  to  adapt  for  H-vertical  polarization.  In 
the  original  design,  the  coax  was  straight. 
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available.)  The  balanced  two-wire  output  from  the  sensor  was  measured  using  a 
digital  ballun  whereby  each  side  of  the  voltage  was  sequentially  recorded  with 
the  other  terminated  in  50  a  and  the  resultant  obtained  by  vector  subtraction. 

The  MGL-9(R)  probe  was  applicable  to  measure  the  vertical  H-field  but 
could  not  be  used  for  the  horizontal  field  component  measurements.  There  we 
used  our  own  make  which,  due  to  its  0.020  inch  semirigid  cable,  could  be  bent 
to  measure  the  horizontal  component  or  the  vertical  component  needed  for  the 
calibration. 

2.  MEASUREMENTS  AND  DATA 

For  all  the  orientations  of  the  model  that  we  measured,  the  axis  of  the 
sensor  loop  was  always  parallel  to  the  fuselage,  thus  the  (fuselage)  axial 
magnetic  field  component  was  measured.  For  each  excitation,  measurements  were 
made  for  nine  probe  positions  relative  to  the  aircraft  as  specified  in  Figure 
2.  The  measurement  procedures  and  data  analyses  follow  and  are  similar  to 
those  used  in  the  ground  plane  tests  (Sec.  3.3).  A  total  of  six  cases  were 
studied,  four  for  the  case  of  H-parallel  to  the  fuselage  and  two  for  the  E- 
parallel  to  the  fuselage. 

Recall  (Fig.  1)  that  in  the  chamber  HQ  is  vertical  and  E1  is  horizontal 
and  the  probe  is  suspended  from  the  ceiling.  For  the  case  of  H-parallel  to 
the  fuselage  the  model  was  supported  vertically  on  its  tail  on  the  pedestal 
and  moved  through  its  nine  measurement  positions.  The  photograph  in  Figure  13 
shows  the  measurement  at  position  9,  top  incidence  parallel  to  the  fuselage. 
The  calibration  measurement  vQ(fn)  is  also  the  reference  measurement  v^(fn)  in 
this  case,  and  is  made  with  the  model  removed.  Extreme  care  was  taken  not  to 
touch  or  disturb  the  model  in  the  process. 

For  the  case  of  E-parallel  to  the  fuselage  the  model  was  placed  with  the 
fuselage  horizontal  on  the  Styrofoam  column  and  again  repositioned  through  the 
nine  measurement  positions.  For  the  reference  measurement  ( f n )  the  model 
was  taken  away,  but  for  the  calibration  of  the  incident  field  the  probe  had  to 
be  bent  to  the  shape  shown  in  Figure  12  to  measure  the  incident  field 
vQ(fn).  Actually,  the  current  on  a  3-in-dia  sphere  was  measured  and  then 


Figure  13. 


Photograph  of  the  model  set  up  in  the  chamber  to  measure  the 
field  at  Pos.  9.  A  paper  scale  is  attached  to  the  model  to 
locate  the  measurement  positions.  The  probe  is  MGL-9(R); 
illumination  Top,  H-parallel  to  fuselage. 
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corrected  for  the  sphere  response  to  obtain  the  incident  field  HQ,  a  technique 
developed  and  used  in  free  space  surface  field  measurements  (Ref.  6). 


Table  2  below  shows  the  illumination  situations  measured  and  the  data 
presented.  The  six  data  sets  follow.  Each  data  set  contains  information  on 
the  excitation,  data  files  used,  errors  E}  and  E2,  and  the  response  curves  for 
sensor  positions  4  and  6.  Again,  in  processing  the  data  the  phase  K  was 
adjusted  for  min  E^  (Sec.  3.3)  to  compensate  for  the  phase  drift  in  the 
instrumentation. 


TABLE  2.  SUMMARY  OF  FREE  SPACE  MEASUREMENTS 
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V.  INTERPRETATION  OF  RESULTS 


i 

I 

I 

'  In  Sections  3  and  4  the  errors  caused  by  the  field  scattered  from  the 

aircraft  were  deduced  for  each  of  the  nine  points  measured  in  the  frontal  area 
of  the  aircraft  for  17  different  excitations  that  are  described  by  direction 
of  incidence  and  polarization  with  respect  to  the  aircraft.  Of  these,  11 
excitation  cases  came  from  the  ground  plane  measurement  study  (Sec.  3)  and  six 
from  the  free  space  measurement  study  (Sec.  4).  In  the  case  of  the  ground 

plane  measurements  the  situations  simulated  were  Nose-on  (e  =  180  deg  ,  Fig. 
11),  Top  (e  3  90  deg),  Tail-n  (e  =  0  deg),  below  Tail-on  (9  =  -10  deg),  and 
above  Tail-on  (e  =  +10  deg).  The  polarization  was  always  E-perpendicular  to 
the  fuselage.  In  addition,  for  the  nose-on  and  top  incidences,  the  measure¬ 
ments  were  made  with  the  ground  plane  rotated  to  20,  45,  and  60  deg.  Such 
rotation  simulates  a  situation  that  is  equivalent  to  two  plane-waves  impinging 
symmetrically  on  the  opposite  sides  of  the  aircraft. 

In  the  free  space  measurement  study,  the  measurements  were  made  for  top, 
bottom,  side,  and  below  side  incidence  with  E-parallel  and  E-perpendicular  to 
the  fuselage. 

The  errors  caused  by  the  scattering  of  the  electromagnetic  field  from  the 
aircraft  were  determined  by  computing  the  average  (Eq.  1)  and  the  RMS  (Eq.  2) 
deviations  between  two  signals:  one  measured  using  a  loop  probe  with  the 
model  present  (test  measurement)  and  the  other  with  the  model  removed 
(reference  measurement). 

Scanning  through  the  computed  error  tables  as  well  as  the  data  curves 
themselves,  one  can  see  that  the  F-106  aircraft  produces  little  scattering  in 
the  forward  region.  The  errors  follow  a  pattern  one  might  expect,  the  largest 
ones  being  observed  near  the  bulkhead  and  decreasing  as  one  goes  away,  or 
further  forward  of  the  bulkhead.  The  errors  also  depend  on  the  direction  of 
arrival  and  the  polarization  of  the  incident  signal  relative  to  the  model. 
They  are  lowest  for  nose-on  and  tail -on  incidences  (E^  =  1  percent)  and 
largest  for  side  incidences  (E^  s  3  -  7  percent).  The  largest  errors  were 
observed  for  the  case  E-vertical  Side-on  incidence  at  -30  deg  below 
horizontal.  This  is  one  of  the  cases  that  simulates  situations  of  the  VPf) 
II/F-106  fly-by  tests  performed  at  Kirtland  Air  Force  Base  in  February  of 
1984. 
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Since  the  errors  measured  were  small  and  approach  the  system  noise  in 
many  cases,  we  have  found  it  difficult  to  deduce  the  optimum  location  for  the 
sensor.  The  system  noise  can  be  attributed,  in  part,  to  physical  changes  in 
the  chamber  dimensions,  possible  bumping  of  the  probe  by  the  operator  when 
repositioning  the  model,  or  changes  in  gain,  frequency,  etc.,  in  the  elec¬ 
tronics  during  the  measurement  sequence.  An  actual  study  to  determine  these 
errors  was  not  performed,  but  we  estimate  that  they  could  affect  the  results 
by  as  much  as  +0.25  percent. 

To  determine  which  of  the  nine  positions  measured  have  the  least  errors 
we  used  two  approaches.  In  the  first  one  we  simply  counted  the  number  of 
times  the  lowest  three  errors  appeared  at  a  given  position  for  all  the  17 
measurements  performed.  Such  a  count  is  shown  in  Table  3.  As  can  be  seen, 
the  lowest  error  appears  seven  times  at  position  9,  four  times  at  position  6, 
and  three  times  at  position  3,  but  because  of  the  (estimated)  0.25  percent 
effect  of  the  system  noise,  the  ranking  given  in  Table  3  is  not  necessarily 
meaningful . 

In  the  next  approach  we  looked  for  the  minimum  error  sensor  position  by 
determining  at  each  position  an  average  error  obtained  by  averaging  over  all 
17  measurement  situations.  Positions  1,  4,  and  7  which  are  closest  to  the 
bulkhead  were  not  included  because  (1)  a  sensor  would  not  normally  be  located 
there,  and  (2)  the  errors  measured  there  are  normally  much  higher  and  may 
erroneously  bias  the  results.  Table  4  gives  a  matrix  of  E1  errors  obtained 
from  the  data  sets  of  Sections  3  and  4.  The  columns  give  errors  at  a  given 
position  for  each  of  the  17  situations.  At  the  bottom  of  each  column  are  the 
average  errors,  and  here  we  see  that  the  lowest  error  (2.34  percent)  is  at 
position  9,  second  lowest  (2.35  percent)  at  position  3,  and  the  third  lowest 
(2.4  percent)  at  position  6. 

Because  the  error  levels  strongly  depend  on  the  direction  and  polariza¬ 
tion  of  illuminating  signals  (see  the  right-hand  column  in  Table  4  that  gives 
the  errors  obtained  by  averaging  over  the  six  measurement  positions),  we  then 
proceeded  to  normalize  the  errors  such  that  the  average  for  each  set  is  2.61, 
the  same  as  the  cumulative  average  for  all  positions  and  illumination  given  in 
Table  4.  Table  5  gives  the  normalized  data,  with  new  averages  computed  for 
each  measurement  position.  Now  positions  3  and  6  show  the  minimum  error  (2.36 
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TABLE  3.  MINIMUM  ERROR  RANKING  FROM  THE  RAW  DATA 


TABLE  4.  AVERAGES  FOR  ACTUAL  DATA 


SET 

E2 

E3 

E5 

E6 

E8 

E9 

EAVG 

3.1 

1.28 

1.44 

1.03 

1.00 

1.00 

1.09 

1.14 

3.2 

2.81 

2.10 

2.88 

1.61 
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0.94 
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1.25 
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4.02 

3.79 

3.61 

3.58 

3.86 

3.5 

3.51 

3.67 

4.05 

4.00 

4.47 

4.22 

3.99 

3.6 

1.42 

1.14 

1.91 

1.25 

1.35 

1.11 

1.36 

3.7 

0.97 

0.87 

0.92 

0.66 

1.38 

0.80 

0.93 

3.8 

0.82 

0.64 

0.82 

0.79 
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0.86 

0.83 

3.9 

1.98 

1.62 

1.68 

1.40 

2.38 
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1.79 
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1.71 

1.27 

2.86 

2.42 

2.65 

1.90 

2.13 

3.11 

2.60 

2.30 

2.55 

2.27 

2.96 

2.25 

2.49 

4.1 

1.96 

1.57 

2.10 

1.96 

2.45 

2.05 

2.02 

4.2 

1.97 

1.78 

1.92 

1.70 

2.24 

1.85 

1.91 

4.3 

4.43 

3.84 

4.38 

3.86 

4.31 

3.81 

4.11 

4.4 

11.40 

7.12 

7.77 

8.39 

8.47 

6.87 

8.34 

4.5 

4.22 

3.56 

5.16 

2.45 

3.63 

2.80 

3.64 

4.6 

2.72 

2.23 

2.49 

2.39 

2.11 

1.95 

2.31 

AVG 

2.89 

2.35 

2.79 

2.41 

2.85 

2.34 

2.61 
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TABLE  5.  AVERAGES  FOR  NORMALIZED  DATA 


F2 

F3 

F5 

F6 

F8 

F9 

FAVG 

2.92 

3.29 

2.35 

2.29 

2.29 

2.49 

2.61 

3.08 

2.30 

3.16 

1.77 

3.55 

1.77 

2.61 

2.30 

2.52 

2.34 

2.44 

2.98 

3.05 

2.61 

2.93 

2.58 
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2.61 
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2.65 

2.61 

2.92 
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2.61 

2.71 
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2.12 

2.61 

2.71 
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3.85 
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2.61 
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2.71 

2.61 

2.88 
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3.46 
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2.61 
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2.61 
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2.03 

2.71 

2.53 
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2.61 

2.69 

2.43 
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3.06 
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2.01 
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2.51 
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2.69 

2.37 

2.19 

2.61 

2.76 
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2.96 
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2.61 

84 


percent)  and  the  positions  9  is  the  third  lowest  (2.39  percent).  Results  of 
this  analysis  support  the  theoretical  expectation  by  Giri  and  Baum  (Ref. 4) 
which  says  that  the  minimum  errors  should  be  expected  at  position  3. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  study  has  shown  that  the  errors  caused  by  the  scattering  from  the 
F-106  aircraft  are  for  all  practical  purposes  negligible  when  the  B-dot  sensor 
is  placed  2.7  m  or  further  ahead  of  the  bulkhead.  Because  of  smallness  of  the 
aircraft  scattering,  the  presence  of  measurement  noise  caused  by  reflections 
from  the  chamber,  equipment  drifts,  etc.,  made  it  difficult  to  clearly  predict 
the  optimum  location  of  the  sensor.  More  analyses  on  the  data  taken  could  be 
carried  out  using,  for  example,  only  the  excitations  (polarizations  and  direc¬ 
tions  of  incidence)  appropriate  for  specific  lightning  situations. 

The  errors  should  also  be  assessed  in  the  time  domain  for  a  specified 
class  of  pulse  shapes.  Sine  the  measured  frequency  domain  data  are  already 
in  digital  form,  such  analyses  would  not  be  difficult  and  could  be  carried  out 
by  first  multiplying  the  (measured)  transfer  functions  by  the  provided  pulse 
(lightning)  spectrum  and  then  transforming  to  the  time  domain.  A  more  realis¬ 
tic  error  value,  such  as  an  error  in  the  peak  value  of  the  field  measured  by 
the  sensor,  could  thus  be  obtained. 


APPENDIX  A 


SEM  STUDIES  ON  THE  F-106 


1.  INTRODUCTION 

The  singularity  expansion  method  as  applied  to  electromagnetic  scattering 
is  a  method  of  characterizing  the  response  of  a  target  in  terms  of  transfer 
function  parameters  such  as  poles  and  residues.  To  this  end,  the  poles  are 
assumed  to  be  invariant,  while  the  residues  vary  with  the  particular  measure¬ 
ment  performed  and  the  type  (including  direction  and  polarization)  of  illumi¬ 
nation.  At  least  in  theory,  the  accurate  location  of  a  few  poles  could  assist 
in  determining  the  identity  of  the  target,  but  in  practice  the  supposedly 
invariant  poles  tend  to  move  about  with  changes  in  the  target's  orientation 
and  illumination.  If  the  noise  and  clutter  are  nondetermini stic,  and  if  it  is 
feasible  to  repeat  the  measurement,  both  the  noise  and  clutter  can  be  reduced 
by  repeating  the  measurement  several  times  and  averaging  the  associated  sig¬ 
nals.  If  the  usual  assumptions  about  the  noise  and  clutter  (such  as  distribu¬ 
tions  which  have  zero  mean  and  are  independently  and  identically  distributed) 
are  satisfied,  combining  K  measurements  will  result  in  an  increase  in  the 
si gnal-to-noise  ratio  by  a  factor  of  K.  By  taking  K  sufficiently  large,  the 
data  may  be  quieted  enough  to  guarantee  accurate  pole  and  residue  extrac¬ 
tion.  Of  course,  in  practice,  obtaining  a  large  number,  K,  of  measurements, 
may  be  somewhat  expensive,  or  in  some  situations,  impossible.  For  example,  if 
measurements  are  made  of  a  moving  target  such  as  an  airplane,  the  requirement 
that  the  measurement  be  repeatable  can  be  satisfied  approximately  at  best,  and 
then  only  for  a  limited  amount  of  time.  Conversely,  if  the  measurement  of  a 
stationary  object  is  repeated,  any  stray  reflections  from  other  stationary 
objects  (or,  in  the  case  of  near-field  measurements,  probe-target  interaction) 
will  be  repeated  exactly.  Under  these  circumstances,  the  clutter  is  determin¬ 
istic  and  the  energy  associated  with  the  clutter  cannot  be  reduced  by  repeat¬ 
ing  and  averaging  the  measurements.  On  the  other  hand,  if  the  noise  and 
clutter  are  nondetermini  Stic,  the  signal-to-noise  ratio  will  increase  as 
additional  data  sets  are  incorporated  into  the  composite  data  set.  It  is 
therefore  desirable  to  include  not  only  similar  measurements,  which  may  be 
combined  using  strict  waveform  averaging,  but  also  those  associated  with 
different  illuminations  of  the  target. 


This  Appendix  presents  an  extension  of  measurement  averaging  that  elimi¬ 
nates  the  requirement  of  identical  measurements  and  thus  permits  further 
increases  in  signal-to-noise  ratio  by  increasing  the  number  of  data  sets  which 
may  be  properly  combined.  A  discussion  of  the  accuracy  of  this  method  is 
given  in  Ksienski  (Ref.  Al) 


2.  FORMULATION 


Using  the  standard  pole  residue  expansion,  the  data  as  a  function  of  the 
circular  frequency  u>  may  be  represented  as 
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where  k  is  an  index  to  the  various  measurements  and  ( jto )  represents  noise 
and  clutter.  We  wish  to  combine  K  such  data  sets  in  a  manner  so  as  to  have  a 
maximum  signal-to-noise  ratio  in  the  composite  data  set.  Assuming  a  linear^ 
combination  with  arbitrary  complex  weighting  coefficients,  the  composite  data 
set  is 
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and  by  interchanging  the  order  of  summation  and  noting  the  invariance  of  the 
poles  sm  with  measurement  number  k, 
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In  general,  it  will  only  be  possible  to  maximize  one  bm  at  a  time.  This  is 
because  for  different  choices  of  m  the  akm  will  not  vary  in  unison  with  k. 
Thus,  the  increased  signal-to-noise  ratio,  and  hence  increased  accuracy,  will 
only  be  obtained  for  the  pole  sm.  However,  since  the  specification  of  the 
pole  is  arbitrary,  several  poles  may  be  obtained  through  successively  emphas¬ 
izing  different  bm.  For  the  present,  we  restrict  ourselves  to  obtaining  an 

o 

improved  estimate  of  s^,  which  necessitates  maximizing  |b^j  .  From  (Eq.  A4), 
we  may  write 
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where  w  =  [w^,W2»...  ,wK]T  and  a  =  [a^,  a2^ , . . .  ,aK^]^.  Without  loss  in  gener¬ 
ality,  the  weighting  vector  may  be  constrained  to  w^w  =  1.  Then,  noting  the 
dot  product  formulation  of  (Eq.  A5)  the  maximum  b^  will  occur  for 
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and  for  this  optimal  choice  of  the  weighting  vector,  b1 


To 


evaluate  the  energy  associated  with  the  noise  N(jo>),  we  take  the  expectation 
of  the  noise  squared 
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and  assuming  the  variance  of  the  noise  is  equal  to  a  in  each  of  the  original 
measurements , 
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Thus  the  energy  associated  with  the  noise  will  not  increase  with  w  normalized 
as  in  (Eq.  A6).  If  K  data  sets  are  combined  each  containing  an  approximately 
equal  excitation  of  the  mode  associated  with  the  desired  pole,  the  signal-to- 
noise  ratio  will  increase  by  a  factor  of  K.  For  any  set  of  residues,  the 
resulting  weighting  coefficients  will  produce  the  maximum  possible  increase  in 
si gnal -to-noise  ratio. 


3.  IMPLEMENTATION 

The  formulation  assumes  initial  estimates  of  residues  in  order  to  compute 
the  weighting  coefficients.  The  algorithm  is  fairly  insensitive  to  errors  in 
these  estimates  of  the  residues  and,  in  fact,  a  single  complex  error  factor 
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associated  with  all  of  the  residues  will  not  adversely  affect  the  weighting 
vector.  The  effect  of  variations  in  the  real  part  of  the  pole  is  minimized  by 
normalizing  the  residues  by  the  real  part  of  the  associated  poles,  i.e.,  (Eq. 
6)  is  used  with  a  redefined  as 


where  s^  refers  to  the  estimate  of  the  pole  s^  obtained  from  the  ith  data 
set.  After  the  composite  data  set  is  formed,  it  must  be  subjected  to  a  pole 
and  residue  extraction  procedure.  The  algorithm  due  to  Levy  (Ref.  A2)  and 
Sanathanan  and  Koerner  (Ref.  A3)  assumes  that  the  poles  and  residues  exhibit 
conjugate  symmetry.  This  is  generally  appropriate  as  it  is  equivalent  to 
requiring  the  frequency  data  to  correspond  to  a  real  function  of  time.  The 
Fcomp  does  not  have  conjugate  symmetry.  Recently,  Tao  and  Zunde  (Ref. 
A4)  extended  the  method  of  Levy  (Ref.  A2)  to  distributions  of  poles  and  zeroes 
which  have  no  symmetry  constraints.  Although  the  algorithm  of  Tao  and  Zunde 
does  accommodate  the  data  set,  it  is  somewhat  more  general  than  required. 
Specifically,  it  does  not  require  the  poles  to  occur  in  conjugate  pairs,  but 
as  may  be  seen  from  (Eq.  A 2),  Fcomp  ( jtu)  does  preserve  the  conjugate  symmetry 
of  the  pole  pairs.  To  more  closely  match  the  known  parametrization  of  the 
data,  an  algorithm  was  produced  which  relaxes  the  conjugate  symmetry  con¬ 
straint  on  the  residues  while  maintaining  the  conjugate  symmetry  constraint  on 
the  poles,  and  the  derivation  is  given  in  Ksienski  (Ref.  Al). 

The  implementation  of  the  algorithm  requires  negative  frequency  informa¬ 
tion  which  is  obtained  from  the  original  data  sets  by  employing  conjugate 
symmetry.  After  the  accuracy  of  the  pole  location  is  improved,  the  next  step 

is  to  refer  the  residues  to  the  new  pole  estimate.  Constraining  a  pole  to  a 

particular  location  is  not  an  easy  task  (Pond  and  Senior,  Ref.  A5),  but  the 
present  algorithm  provides  a  simple  alternative.  Since  an  original  data  set 
corresponds  to  a  real  function  of  time,  it  is  an  even  function  of  frequency. 
A  second  set  of  data  which  is  an  odd  function  of  frequency  may  be  added  to  the 
original  data  and  the  sum  expanded  as  a  series  of  poles.  Noting  that  the 

poles  are  constrained  to  occur  in  conjugate  pairs,  the  component  of  the  summa¬ 
tion  associated  with  the  original  data  may  be  obtained  by  retaining  the  even 
part  of  the  residues.  Both  the  original  data  and  the  second  set  are 


constrained  to  have  the  same  poles  in  the  expansion  and,  by  appropriate 
placement  of  the  poles  associated  with  the  second  set  of  data,  it  is  possible 
to  constrain  the  pole  locations  in  the  expansion  of  the  original  data.  The 
residues  are  then  obtained  relative  to  these  new  constrained  pole  locations 
without  any  apparent  adverse  side  effects.  An  additional  advantage  is  that 
the  accuracy  of  the  location  of  the  adjacent  pole  is  also  improved. 

4.  RESULTS 

It  was  desired  to  extract  the  two  poles  associated  with  the  fuselage  and 
wing  resonances.  Casual  inspection  of  the  data  indicated  that  these  poles 
would  have  imaginary  parts  of  approximately  6  and  14  MHz,  respectively.  The 
algorithm  described  above  was  applied  to  14  data  sets  (Liepa  and  Pennock,  Ref. 
A6)  which  were  selected  as  having  relatively  large  components  associated  with 
the  desired  poles.  The  data  sets  used  were  FT1476,  FT1517,  FT1601,  FT1840, 
FT1919,  FT1809,  FT1801,  FT1765,  FT2465,  FT2433,  FT2372,  FT1757,  FT1617,  and 
FT1309.  Figures  A1  through  A4  show  plots  of  these  data.  The  poles  and  resi¬ 
dues  which  were  extracted  are  listed  in  Table  Al. 

From  Reference  A. 6  are  presented  here  Tables  A2,  A3  and  Figure  A5  that 
identify  each  data  set  with  measurement  situations. 
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TABLE  Al.  POLES  AND  RESIOUES  FOR  NASA  F-106  AIRCRAFT 


^N>v>Poles 

Residue 

-1.272+j6.457 

-1 . 459+j 1 4.45 

- 2. 609+ j 28. 66 

FT1465 

-1. 291 -j. 2784 

4.153-j2.302 

1 . 665- j 5. 149 

FT1309 

2562+ j. 07 608 

1 .805-j .6113 

-.223 5+j. 09827 

FT! 517 

1 430+ j. 1857 

-. 4043+ j. 5228 

-. 4470- j. 3472 

FT1601 

-1.163+j4.238 

. 7896- j. 5094 

-.  3829+ j. 2955 

FT1617 

-21166+j3.795 

-3. 656- j2. 854 

-2. 422- j. 11 26 

FT!  7  57 

-9 - 230+j 18. 22 

1 .214-jl .829 

-9.374+j4.474 

FT! 765 

-3.907-j8.879 

-. 2404+ j. 02450 

-1 .71 6— j7 . 743 

FT1801 

-6.985+J2.436 

2.552+jl .544 

-3. 542+ j 9. 433 

FT1809 

4.911-j20.30 

. 5698+ j. 8693 

1.887-j4.021 

FT1849 

-.1 381 +j. 2256 

-.04788+j. 09923 

. 4091 + j 1.1 93 

FT1917 

- 2. 084- j 1.449 

.  1 096+j . 31 68 

1. 280+ j. 4045 

FT2373 

-. 01 067+j. 07738 

- . 03776- j . 1 230 

1.074-j2.020 

FT2433 

3.453+jl .439 

.861 5+j .21 69 

2. 048- j 1.450 

FT2465 

3. 126+jl .938 

.6381-j .2114 

2. 608+ j 2. 124 

m  tv  m  uw 


TABLE  A2.  MEASUREMENT  MATRIX 
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APPENDIX  B 

LIGHTNING  SIMULATION  STUDIES 

Recent  measurements  of  the  electric  and  magnetic  fields  of  lightning 
processes  reveal  that  the  fields  can  vary  on  a  submicrosecond  time  scale.  The 
rapid  variation  of  the  electromagnetic  fields  implies  that  the  currents 
involved  in  a  lightning  discharge  also  can  vary  on  a  submicrosecond  scale 
(Refs.  Bl,  B4,  B3).  Airborne  measurements  have  directly  confirmed  that  cloud 
discharge  processes  can  have  submicrosecond  current  rise  times  (Refs.  Bl,  B5, 
B2) . 

Lightning  processes  with  submicrosecond  current  and  field  rise  times  are 
a  hazard  to  aircraft  because  aircraft  resonances  can  be  efficiently  excited. 
Direct  lightning  strikes,  as  well  as  near  misses,  can  excite  resonances  which 
pose  a  threat,  especially  to  newer  aircraft.  The  increased  use  of  nonmetal  lie 
structural  materials  has  reduced  the  electromagnetic  shielding  of  the  interior 
of  the  aircraft,  making  flight-critical  low-voltage  digital  electronics  sus¬ 
ceptible  to  interference  or  damage  (Ref.  B6). 

An  aircraft  in  flight  can  become  attached  to  a  lightning  channel  by 
initiating  discharges  which  connect  to  the  stepped  leader  or  to  any  other 
phase  of  the  discharge  in  which  there  is  a  propagating  channel.  In  addition, 
an  aircraft  can  cause  a  lightning  discharge  which  would  not  have  occurred  if 
the  aircraft  had  not  been  there.  The  lightning  discharge  is  triggered  by 
field  enhancement  caused  by  the  presence  of  the  aircraft  in  a  region  where  the 
electromagnetic  field  is  already  high,  resulting  in  the  generation  and  propa¬ 
gation  of  leaders  in  two  directions  away  from  the  aircraft  toward  regions  of 
opposite  charge.  Charge  on  the  aircraft  due  to  precipitation  interactions  can 
increase  the  probability  of  a  lightning  strike.  Charge  deposited  in  the  wake 
of  the  aircraft  due  to  precipitation  charging  can  guide  the  lightning  to  the 
plane  or  shield  the  plane  from  it.  An  extensive  discussion  of  the  lightning- 
aircraft  interaction  is  given  in  the  paper  by  Clifford  and  Kasemir  (Ref.  B7). 

The  attachment  of  a  typical  lightning  discharge  to  an  aircraft  in  flight, 

as  described  in  the  paper  by  Uman  and  Krider,  happen  as  follows:  "A  lightning 

_  1. 

leader  with  a  typical  charge  density  of  10  c/m  approaches  the  aircraft  at  a 

c  c 

velocity  of  10  to  10  m/s.  When  the  fields  on  the  wings,  nose,  tail,  and 
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other  extremities  exceed  a  critical  value,  corona  will  begin  and  outward-going 
discharges  will  be  initiated.  These  discharges  will  lower  the  electric  field 
at  the  surfaces  of  the  plane  because  of  shielding.  The  'striking  distance'  is 
determined  by  the  distance  to  the  incoming  leader  at  the  time  when  the  air¬ 
craft  initiates  the  outward-going  discharges.  The  time  for  the  leader  and 
connecting  dicharge  to  join  is  probably  of  the  order  to  10  ys,  and  the  strik¬ 
ing  distance  is  probably  of  the  order  of  a  plane  length.  When  contact  takes 

Q  Q 

place,  the  aircraft  will  be  raised  to  the  10  to  10*  V  potential  of  the  light¬ 
ning  channel  in  a  characteristic  charging  time  determined  by  the  channel  surge 

impedance,  roughly  1000  a,  and  the  aircraft  capacitance  10“9  to  10"10  °F. 

7 

Thus  the  charging  time  will  be  10  to  10“°  s,  and  the  rate  of  change  of  the 

o  IQ 

aircraft  voltage  will  be  10°  10  V/ys.  Contact  will  also  cause  the  electri¬ 
cal  field  at  the  surface  of  the  plane  to  reverse  direction,  and  the  aircraft 
may  produce  additional  corona  and  leaders.  In  any  event,  the  leader  current 
will  propagate  through  the  aircraft  in  a  relatively  steady  fashion,  and  then 
there  may  be  superimposed  fast  current  pulses  due  to  leader  steps,  K-changes, 
return  strokes,  or  any  other  implusive  lightning  process."  (Ref.  86). 

The  purpose  of  this  investigation  was  to  investigate  methods  of  analyzing 
and  simulating  in  the  laboratory  a  direct  lightning  strike  on  an  aircraft  in 
the  frequency  domain.  The  description  of  the  lightning  attachment  process  and 
the  presence  of  submicrosecond  components  in  lightning  currents  and  fields 
suggest  that  the  primary  effect  of  a  lightning  attachment  to  an  aircraft  in 
the  frequency  domain  would  be  to  excite  and  modify  natural  aircraft  resonances 
and  possibly  produce  new  resonances  associated  with  the  attachment  of  the 
highly  conducting  lightning  channel. 

The  natural  resonant  frequencies  of  an  aircraft  can  be  evaluated  experi¬ 
mentally  by  illuminating  a  scale-model  of  the  aircraft  in  an  anechoic  chamber 
with  a  plane  wave  of  appropriately  scaled  frequency  and  measuring  the  surface 
currents  i  duced  on  the  model.  We  felt  that  the  effect  of  a  lightning  attach¬ 
ment  on  the  natural  resonant  frequencies  of  an  aircraft  could  be  evaluated 
experimentally  by  attaching  a  structure  simulating  the  lightning  channel  to 
the  scale  model  of  the  aircraft,  placing  the  model  in  the  chamber  and  measur¬ 
ing  the  surface  currents  induced  on  the  surface  of  the  model  airplane  when  it 
is  illuminated  by  a  plane  wave. 
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Our  primary  task  thus  lay  in  choosing  an  appropriate  lightning  model  and 
finding  a  structure  (a  wire,  lossy  string,  etc.)  which  would  accurately  repre¬ 
sent  the  lightning  model  in  the  measurement  frequency  range.  The  structure 
attached  to  the  aircraft  scale-model  could  then  be  used  to  investigate  the 
effect  of  lightning  attachment. 

A  model  can  be  defined  as  a  physical  or  mathematical  construct  which 
approximates  to  some  degree  certain  aspects  of  natural  or  man-made 
phenomena.  There  are  basically  three  levels  of  sophistication  in  the  mathema¬ 
tical  modeling  of  lightning. 

The  simplest  lightning  channel  model,  which  provides  the  closest  approxi¬ 
mation  to  nature,  assumes  a  temporal  and  spatial  form  for  the  channel  current 
and  then  uses  the  current  to  calculate  the  remote  fields.  A  more  sophisti¬ 
cated  approach  involves  mathematically  describing  the  lightning  channel  as  an 
RLC  transmission  line  with  distributed  circuit  elements  that  may  vary  with 
height  and  time.  The  intent  there  is  to  predict  a  channel  current  as  a  func¬ 
tion  of  height  and  time,  and  to  use  this  current  to  calculate  the  fields.  The 
most  sophisticated  modeling  method  attempts  to  describe  the  detailed  physics 
of  the  lightning  channel  using  equations  of  conservation  of  mass,  momentum, 
and  energy,  equations  of  state  and  Maxwells'  equations  (Ref.  B6).  Our 
interest  in  the  several  models  is  to  find  a  valid  description  of  the  lightning 
channel  geometry  which  could  serve  as  the  basis  from  which  to  find  a  reason¬ 
able  experimental  model. 

For  our  study  we  chose  to  expand  on  the  models  used  by  Yang  and  Lee  (Ref. 
38).  These  lightning  models  were  relatively  sophisticated  and  used  channel 
geometries  that  were  amendable  to  experimental  implementation.  In  addition, 
since  the  natural  frequencies  of  a  post  attached  to  a  lightning  return  strike 
were  calculated  analytically  in  the  above  report,  the  natural  frequencies  of  a 
post  attached  to  a  lightning  return  stroke  determined  experimentally  could  be 
compared,  thus  providing  a  method  to  check  the  validity  of  the  experimental 
model . 

Yang  and  Lee  used  three  lightning  models,  of  which  two  were  special  cases 
of  the  third.  The  first  model  considered  was  the  resistive  model  of  Volland 
(Ref.  39)  who  showed  that  a  lightning  channel  and  its  surrounding  corona  could 
be  represented  by  a  finitely  conducting  rod  with  an  effective  radius.  Such  a 
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representation  can  be  considered  a  transmission  line  (Ref.  BIO).  The  required 
parameters  are  the  effective  radius  re  and  the  effective  resistance  per  unit 
length  R'.  In  term?  of  these  two  parameters  the  lightning  channel  has  a 
characteristic  impedance 
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where  s  is  the  complex  frequency,  and  ZQ,  c,  uQ  are  the  free-space  wave  impe¬ 
dance,  speed  of  light,  and  permeability. 

The  corona-sheath  model  has  been  suggested  in  Reference  Bll  for  a  lightn¬ 
ing  channel  and  its  surrounding  corona.  The  lightning  current  flows  only  in 
the  perfectly  conducting  center  channel,  while  all  the  electric  charges  reside 
on  an  effective  corona  surface.  The  required  parameters  for  the  corona-sheath 
model  are  the  effective  radius  of  the  lightning  center  channel  rQ  and  the 
effective  radius  of  the  corona  surface  r£.  The  effective  corona  radius  rc  is 
related  to  the  charge  per  unit  length  of  the  corona  via 


where  Eb  and  eQ  are  the  air  breakdown  electric  field  and  permittivity  of  the 
surrounding  air  respectively.  This  model  can  also  be  considered  as  a 
transmission  line  with  characteristic  impedance 


r  Z  sr  sr  1/2 

z‘(s)  ■£♦*<«),  *c(x)  =  [tn  (-j£)  w 


A  composite  of  the  two  above  models,  the  resisti ve-corona  sheath  model,  is  the 
same  as  the  corona  sheath  model  except  now  the  center  channel  can  have  a 
finite  resistivity. 

The  geometries  of  the  theoretical  lightning  models  described  were 
conducive  to  constructing  experimental  models,  but  materials  from  which  to 
make  the  models  were  not  available.  The  resistive-model  could  be  simulated  by 
a  wire  or  string  of  appropriate  size  and  resistivity.  Wires  were  available 
with  the  properly  scaled  physical  dimensions  but  wires  with  appropriate 


conductivity  were  not.  Yang  and  Lee  evaluated  the  natural  frequencies  of  a 
post  attached  to  lightning  channels  of  several  different  values  of  DC 
resistance  per  unit  length.  The  appropriate  value  of  natural  lightning 
channel  resistance  is  still  a  subject  of  much  research.  Thus,  we  would  like 
to  have  experimental  models  with  several  values  of  channel  resistance. 

When  the  dimensions  and  resistivities  of  the  lightning  channels  evaluated 
in  the  Yang  and  Lee  work  were  appropriately  scaled  to  correspond  to  the  fre¬ 
quency  range  of  our  experimental  facility,  it  was  found  that  wires  with  DC 
resistances  of  0.576  to  5760  a/ft  would  be  necessary  to  duplicate  the  situa¬ 
tions.  The  only  available  wires  of  the  appropriate  diameter  had  DC  resis¬ 
tances  varying  from  ~0  to  16.4  a/ft.  Thus,  our  intent  to  simulate  the 
frequency  domain  effect  of  lightning  using  a  wire  attached  to  a  post  protrud¬ 
ing  from  a  ground  plane  inside  the  chamber  was  impossible  to  achieve. 

Measurements  were,  however,  made  with  the  available  wires  to  see  if  the 
small  variation  in  DC  resistance  had  any  effect  on  the  response  of  the  wires 
in  the  anechoic  chamber.  The  wires  are  listed  in  Table  Bl.  The  same  plane 
(c.f.  Fig.  8,  Sec.  3)  used  in  the  ground  plane  study  of  the  F-106  was 
employed.  One  end  of  each  wire  sample  was  soldered  to  the  center  of  a  coaxial 
protruding  from  the  center  of  the  ground  plane  and  the  other  was  attached  to  a 
metal  rod  which  was  then  stuck  inside  an  absorber  cone  on  the  side  of  the 
chamber.  The  wires  were  illuminated  with  a  plane  wave  whose  electric  vector 

was  along  the  wire  axis.  The  plane  wave  frequency  was  stepped  from  118  to 

4400  MHz  and  the  magnitude  and  phase  of  the  voltage  feeding  the  network 
analyzer  (50  a  system)  was  recorded.  The  wire  was  then  replaced  by  an  MGL-8 
B-dot  sensor  and  the  measurements  repeated.  Dividing  the  wire  measurement  by 
the  MGL-8  measurement  and  making  the  response  correction  for  the  MGL-8  gives 
the  wire  response  in  terms  of  V/HQ,  where  HQ  is  the  incident  magnetic  field. 

A  plot  of  the  response  of  wire  G  with  a  DC  resistance  per  unit  length 

of  ~0  a/ft  is  shown  in  Figure  Bl.  A  plot  of  the  response  of  the  wire  D  with  a 

DC  resistance  per  unit  length  of  16.4  a/ft  is  shown  in  Figure  B2.  From  the 

plots  it  is  evident  that  changing  the  resistance  per  unit  length 

from  ~0  a/ft  to  16.4  a/ft  has  little  effect  on  the  response  of  the  wire.  The 

data  plotted  in  Figures  Bl  and  B2  have  been  digitally  filtered  with  a  routine 

which  preserves  100  cycles  per  plot  to  reduce  noise  in  the  measured  data. 


Examination  of  Figures  81  and  82  shows  that  the  phase  of  the  wire  response 
increases  linearly  with  frequency.  The  amplitude  of  the  wire  response 
decreases  linearly  from  118  MHz  to  approximately  2500  MHz  and  then  becomes 
fairly  constant.  A  periodic  standing  wave  pattern  is  super-imposed  on  the 
wire  response  and  can  be  related  to  the  wire  length  and  current  reflections  at 
each  end  of  the  wire  even  though  the  tip  of  the  wire  was  embedded  in  the 
absorber  of  the  chamber  wall.  The  responses  of  wires  A  through  G  were  virtu¬ 
ally  the  same. 

As  a  result  of  these  measurements  we  could  essentially  only  implement  the 
resistive  model  in  the  case  of  a  conducting  channel  of  infinite 
conductivity.  To  go  further  it  is  necessary  to  find  new  wire  materials  or 
other  structures  with  which  to  implement  the  corona-sheath  and  resistive 
corona-sheath  models.  It  has  been  suggested*  that  plastic  tubing  filled  with 
a  conducting  electrolyte  solution  could  be  suitable  for  the  corona-sheath 
model  since  the  channel  resistance  could  be  changed  by  changing  the 
electrolyte  solution. 

Evaluating  the  effect  of  a  lightning  attachment  on  the  natural 
frequencies  of  an  aircraft  in  the  frequency  domain  experimentally  using  scale 
models  in  the  anechoic  chamber  shows  promise.  However,  more  work  needs  to  be 
done  to  find  valid  experimental  models  of  lightning  channels. 


TABLE  81.  LIST  OF  WIRES  TESTED 


Wire  Composition/ All oj 


Ni chrome 


Gauge  Diameter  (mm) 


0.4064 


Nichrome 


0.2286 


10.58 


Chromel -A 


0.3302 


Chromel -A 


0.1778 


A1 umel 


0.3302 


Kanthal 


0.3302 


Copper 


0.7112 


*  Baum,  C.E.,  personal  communication.  May  1984. 
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Figure  81.  Voltage  received  by  22  gauge  copper  wire  (*0  n/ft) 
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Ffgure  82. 


Voltage  received  by  24  gauge  chromel-A  wire  (-  16.4  a/ft). 
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